DISCONTINUITIES OF THE INTEGRATED DENSITY OF 
STATES FOR RANDOM OPERATORS ON DELONE SETS 

STEFFEN KLASSERT *■•*, DANIEL LENZ 2 -* AND PETER STOLLMANN 3 <* 



1 Fakultat fur Mathematik, Technische Universitat Chemnitz, D-09107 Chemnitz, 
Germany; E-mail: S.Klassert@mathematik.tu-chemnitz.de, 



2 Fakultat fur Mathematik, Technische Universitat Chemnitz, D-09107 Chemnitz, 
Germany; E-mail: D.Lenz@mathematik.tu-chemnitz.de, 



3 Fakultat fur Mathematik, Technische Universitat Chemnitz, D-09107 Chemnitz, 
Germany; E-mail: P.Stollmann@mathematik.tu-chemnitz.de 



Abstract. Despite all the analogies with "usual random" models, tight bind- 
ing operators for quasicrystals exhibit a feature which clearly distinguishes 
them from the former: the integrated density of states may be discontinuous. 
This phenomenon is identified as a local effect, due to occurrence of eigenfunc- 
tions with bounded support. 



1. Introduction 



In the present article we study the occurrence of discontinuities in the integrated 
density of states (IDS). For the special case of a tight binding model associated with 
the Penrose tiling the occurrence of this effect has been known since quite some time 
as witnessed for example in |ATF, FATK, KF, KS|. We present two results. The 
first aims at showing that the occurrence of jumps in the IDS cannot be excluded 
by global assumptions concerning, e.g. ergodic or combinatorial properties. To 
this end we present a theorem saying that starting from some model of aperiodic 
order (phrased in the language of Delone sets) one can construct a model that is 
" basically the same" and gives rise to a tight binding operator with a discontinuity 
in the IDS. Here "basically the same" is cast in the notion of "mutually locally 
derivable" for Delone dynamical systems. We discuss this notion analogous to the 
respective notion for tilings found in BSJ]. In the construction we use that the 
laplacians on certain graphs have finitely supported eigenfunctions. It now becomes 
clear that it is the more complex structure of graphs in higher dimension that makes 
such a phenomenon possible. On lattices (and, consequently, in one-dimensional 
systems) such finitely supported eigenfunctions cannot occur. 

Our second theorem says that this is the only possibility to create a jump of 
the IDS, at least when starting from a reasonable Del one d ynamical system. It is 
a consequence of a rather strong ergodic theorem in LS3 ]. This theorem states 
that the IDS is in fact the uniform limit of eigenvalue counting distributions. It is 
therefore substantially stronger than the weak convergence results typically proven 
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in connection with the IDS see for example [ BLT| , ||, Jk|. The fact that such a 



strong convergence holds true is special for models of aperiodic order and not met 
in usual random systems. 

2. Notation and results 
In this section we introduce some notation and present our results. We will use 



the same setting as the one found LS2|. For completeness reasons we include the 
necessary definitions. 

Let d > 1 a fixed integer and all Delone sets, patterns etc. will be subsets of M. d . 
The Euclidean norm on W 1 will be denoted by || • ||. For r £ 1R + and p G M. d , we let 
B(p,r) be the closed ball in W l centered at p with radius r. 

A subset ui of M. d is called Delone set if there exist r(ui) and R(u>) > such that 
2r(cj) < \\x — y\\ whenever x,y € lu with x ^ y, and B(x,R(uj)) flw ^ for all 

x g R d . 

We are dealing in this paper with local structures of Delone sets, therefore the 
restrictions of u) to bounded subsets of K d are of particular interest. In order to 
treat these restrictions, we introduce the following definition. 

Definition 2.1. (a) A pair (A, Q) consisting of a bounded subset Q of R. d and 
AcQ finite is called pattern. The set Q is called the support of the pattern, 
(b) A pattern (A, Q) is called a ball pattern if Q = B(x 7 r) with x G A for suitable 
x G M. d and r G (0, oo). 

The diameter and the volume of a pattern are defined to be the diameter and 
the volume of its support respectively. 

We will have to identify patterns which are equal up to translation. More 
precisely, on the set of patterns we introduce an equivalence relation by setting 
(Ai,Qi) ~ (A 2 ,Q2) if and only if there exists a t G M. d with Ai = A 2 + t and 
Qi = Q2 + t. The class of a pattern (A, Q) is then denoted by [(A, Q)]. Obviously 
the notions of diameter, volume occurrence etc. can easily be carried over from 
patterns to pattern classes. 

Every Delone set w gives rise to a set of pattern classes, P{oj) = {Q A uj : 
Q C R d bounded and measurable}, and to a set of ball pattern classes Vb{oj) = 
{[B(p, r) A uj] : p G w, r G M + }. Here we set Q A uj = (uj n Q, Q). We define the 
radius s = s(P) of an arbitrary ball pattern P to be the radius of the underlying 
ball. For s G (0, 00), we denote by V%(uj) the set of ball patterns with radius s. A 
Delone set is said to be of finite type or of finite local complexity if for every radius 
s > the set V%(uj) is finite. 

The Hausdorff metric on the set of compact subsets of M. d induces the so c alled 
natural topology on the set of closed subsets of M d . It is described in detail in |LS2| 
and shares some nice properties: firstly, the set of all closed subsets of M. d is compact 
in the natural topology. Secondly, and this is of prime importance in view of the 
dynamical system we are to consider, the natural action T of M d on the closed sets 
in R d given by T t G = G + t is continuous. 

Furthermore, a Delone dynamical system (DDS) consists of a set 51 of Delone 
sets, which is invariant under the shift T and closed in the natural topology. A DDS 
is said to be of finite type (DDSF) if Uu 6 nP|(w) is finite for every s and the set 
T{Sl) of patterns classes associated to a DDS VL is defined by P(fl) — U we n'P(w). 
Due to the compactness of the set of all closed sets in the natural topology a DDS 
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17 is compact. We refrain from a precise discussion of the topology but we give the 
following Lemma from JLS2Q . 



Lemma 2.2. If (17, T) is a DDSF then a sequence (tOk) converges to lo in the 
natural topology if and only if there exists a sequence (tk) converging to such that 
for every L > there is an fco £ N with (u>k + tk) fl B(0, L) = lo D B(0, L) for 
k > fco- 

Roughly speaking, w is close to lo if lo equals w on a large ball up to a small 
translation. We now recall some standard notions from the theory of dynamical 
systems and some available equivalent "combinatorial" characterizations. 

A dynamical system (17, T) is called minimal if the orbit {T t LO : t E M. d } of any 
lo is dense in 17. For a DDS this is equivalent to the property that P(£l) = V(lo) 
for any lo. This latter property is called local isomorphism property in the tiling 
framework; see poll . A sequence (Qk) of subsets in W 1 is called a van Hove 



sequence if the sequence \d R Qk\\Qk\~ 1 tends to zero for every R £ (0, oo). Here, 
d R Q denotes the set of those x E R d whose distance to the boundary of Q is 
less than R. Furthermore, a dynamical system (17, T) is called uniquely ergodic 
if it admits only one T-invariant measure (up to normalization). For a Delone 
dynamical system, this is equivalent to the fact, that for every nonempty pattern 
class P the frequency 

u(P) = lim \Q k \-Hp{Q k Aw), 

k — >oo 

exists uniformly in lo e ^ for every van Hove sequence (Qk)- We call a dynamical 
system (57, T) strictly ergodic if it is minimal and uniquely ergodic. Note that in 
this case the frequency v(P) is positive for every P eP 

Definition 2.3. Let (17, T) be a DDSF. A family (A u ) of bounded operators A u : 
1 2 {lo) — ► ^ 2 ( w ) *s called a random operator of finite range if there exists a constant 
rA with 

• Au(x, y) = whenever \\x — y\\ > rA- 



Au(%,y) only depends on the pattern class (B(x,Ta) U-B(y, Ta)) Acj. 



Note that the above defined operators provide a framework including Laplace 
type operators defined on t 2 {JL d ). Let us mention, that for Hi, u> E D. the opera- 
tors act on different Hilbert spaces 1 2 (lo) and £ 2 (lo) unless lo and lo differ only by 
translation. Thus, to deal with operators on £ 2 (lo) is more complicated than in the 
lattice case. 

The aim of this article is to discuss the phenomenon of discontinuities of the 
integrated density of states of random operators (A u ) on a DDSF (S7,T). This 
might rather come as a surprise in view of what is known for random models as 
well one dimensional quasicrystals. It turns out that this phenomenon occurs if 
and only if there exist locally supported eigenfunctions of (A u ). One can find 



examples of locally supported eigenfunctions on the Penrose tiling in [KS| and 



[ATF |. An eigenfunction / is said to be locally supported if supp/ C K, with 



K a compact set. The phenomenon of locally supported eigenfunctions is by no 
means pathological. Rather from any given DDSF (17, T) we can construct an in 
some sense local equivalent DDSF (17 h ,T) such that a random operator of finite 
range (A 1 ^) defined on (17 b ,T) has locally supported eigenfunctions. More precise: 
(17, T) and (17 h ,T) are mutually locally derivable (MLD). The equivalence concept 
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of mutual local derivability for tilings was discussed in detail in [BSJ|. This will all 
be discussed below. Our first result reads as follows. 

Theorem 1. Let (fi,T) be a arbitrary DDSF. Then there exists a DDFS (fi b ,T) 
and a random operator of finite range {A b w ) such that (O, T) and (tt b ,T) are mutu- 
ally locally derivable and (A^) has locally supported eigenf unctions with eigenvalue 
E for every uj G Sl b . Moreover, {A b w ) can be chosen to be the nearest neighbor 
laplacian of suitable graph. 

Note that for a selfadjoint random operator A and bounded Q C K d the restric- 
tion A u \q defined on £ 2 (Q flw) has finite rank. Therefore, the spectral counting 
function 

n(A LU ,Q)(E) :— #{ eigenvalues of A u \q below E} 
is finite and -^n^A^^Q) is the distribution function of the measure Pq", defined 
by ' " 

(p£» )¥ >) := j^rtrMAJg)) for <p G C b (R). 

On a uniquely ergodic DDSF the measures Pq ^ converge in distribution to a 
measure p A , called the integrated density of states (IDS) for any van Hove se- 



quence Qj. as k — > oo. This is described in [LS2, LS3|. There one can also find an 
interpretation of the IDS as a certain trace on a von Neumann algebra. Now, we 
can state our main theorem. 

Theorem 2. Let (f2, T) be a strictly ergodic DDSF. Let A be a selfajoint random 
operator of finite range. Then E is a point of discontinuity of p A if and only if 
there exists a locally supported eigenfunction of A u to E for one (all) lu G fi. 



Remarks 2.4. (1) It rather straightforward to see that locally supported eigen- 
functions lead to a discontinuity of the IDS. The more interesting part of the 
equivalence is discontinuities only happen in that way. 

(2) The theorem holds also in the tiling setting. Here a single tile of the original 
tiling will be replaced by 9 tiles of a new tiling which is MLD to the originally 
given one. 

(3) As pointed out already the theorem gives rise to a complete characterization 
of the phenomenon of locally supported eigenfunctions in quasicrystal settings 
(i.e. DDSF an tiling settings). 

(4) Let us emphasize that the integrated density of states is continuous in the 
case of almost periodic and random operators on lattices. Due to the more 
complex geometry this does not follow in the quasicrystal framework. 

3. Preliminaries 

In this section we will study the equivalence concept of MLD. Further we are 
going to construct a map, which maps any given DDSF to a DDSF which is MLD 
to the given one and where a random operator exists who has locally supported 
eigenfunctions. As well we state some tools which will be used later on to prove 
our results. 



As mentioned already the equivalence concept of MLD for patterns on tilings was 
discussed in [BSJ|. We give the obvious definition for Delone dynamical systems of 
finite type. 
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Definition 3.1. Let (Q,T) and (fl b ,T) a DDSF. A map D : fl -> fl b is called 
a local derivation map if there exists a radius rp> > such that D(u>) Pi {x} — 
(t + D(lu)) fl {x} holds whenever uj D B(x,rp>) = (t + uj) n B{x,Td)- In this case 
(A b ,T) is called locally derivable from (A,T). Two DDSF (fl,T) and (fl b ,T) are 
mutually locally derivable if (fl,T) is locally derivable from (fl b ,T) and vice versa 
(with a possibly different radius r' D ). 

Note that the map D is local in the sense that D(ui) fl B(x, s) only depends on 
uj (~) B(x,s + 2r D ). 

Proposition 3.2. Let (0, T) be a DDSF, D : fl — > fl b , uj i— > -D(w) a Zoca^ derivation 
map. Then D is continuous with respect to the natural topology. 

Proof. This is immediate as D is local and the topology is local in the sense of 



Lemma 2.2. □ 



What we want to do now is to put a well scaled local structure in a given Delone 
set whenever a certain pattern occurs. Let w be a Delone set and P be a ball 
pattern class with P G V- Then, we define ujp to be the set of all occurrences of P 
in R d , i.e. 

(1) lop = {t e R d : B(t, s(P)) A uj = P}. 

Now, let (fl, T) be a DDSF and r < r(uj) for all w e fl, G be a finite graph with 
Vg the set of vertices of G contained in M. d . Furthermore let diam(G) = and 
Vq C B(0, 25). We use this finite graph to define a local derivation map by setting 
D P y G (uj) =ujU {t + V G :t <E tup} for w G fl and fl 6 := {D P y G (oj) : uj g fl}. Then, 

Dp t v G ■ fl — > fl 6 , w !— > -Dp.VgC^) 

is a local derivation with inverse given by the local derivation map 

H P y G : fl b -» fl, H P y G (uj b ) = {x e uj b : uj b D B(x, §) = w b n B(ar, ^)}. 

Note that (fl,T) is also a local derivation of (fl b ,T). Thus, (fl,T) and (fl b ,T) are 
mutually locally derivable. 

Remarks 3.3. Let fl and fl b be as above. Then 

(1) If (fl, T) is a DDSF then (fl b , T) is a DDSF as well. 

(2) If (fl, T) is a uniquely ergodic DDSF, the same holds for (fl b , T). 

(3) The frequency of G in fl b is the same as the frequency v(P) of P in fl. 

The following two ingredients are essential for the proof of Theorem [|. The 
first one is one of the main results from ESj| . It relies on a strong ergodic type 



theorem proven there (see [ Len ] for a study of uniform ergodic theorems in the one 
dimensional case). 

Theorem 3.4. Let (fl,T) be an aperiodic strictly ergodic DDSF. Let A be a self- 
adjoint operator of finite type and (Qk) be a van Hove sequence in R d . Then, the 
distribution functions of converge uniformly to the distribution function of the 
measure p A and this convergence is uniform in uj G fl. 

The second one is a well known dimension argument from linear algebra which 
we will state for completeness reasons. 

Proposition 3.5. Let H be a finite dimensional Hilbert space, U, V subspaces of 
H with dimC/ > dimV then dimV- 1 D U > 0. 
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4. Proofs 

In this section we prove Theorem [l] and Theorem ^. 

Proof of Theorem To prove the theorem we start with the construction of a 
DDSF where a random operator of finite range {AJ) exists which has locally sup- 
ported cigenfunctions. The starting point is a small graph G/j n = (Vfi n , Efi n ) and 
an eigcnfunction u fi n of the associated nearest neighbor laplacian. For dcfinitcncss 
sake consider Figure 1. The values of Uft n are indicated near the corresponding 
vertices. Here the eigenvalue is E = 0. 




Figure 1. The finite graph Gfi 



It is clear that whatever edges reach out of the four corners in a larger graph 
extending Gfi n , the extension of Ufi n by to the larger vertex set will still constitute 
an eigenfunction of the laplacian on the large graph. It is now easy to implement 
this picture into a given DDSF. In fact, let (fl, T) be a DDSF and P be a ball 
pattern. We use the local derivation map Dpy fin discussed in the last section to 
put in Vfi n from above, scaled properly, whenever P appears. It is obvious by 
the definition of Dpy, in , that this gives rise to a DDSF (f2 b ,T) which is locally 
derivable from (fi, T) and vice versa. 

Obviously we get a random operator A b with locally supported eigenfunctions 
by taking for A^ the nearest neighbor laplacian on the copies of Efi n in u> and 
consistent matrix elements otherwise. 

□ 



Remarks 4.1. (1) The simplest case of the construction made above is of course 
given by choosing P = ({x}, B(x, r)) with r < r(u>). Then the graph Gfi n is 
glued at any point of the underlying Delone set. The corresponding £ 2 space 
is just a direct sum (or tensor product) and that applies to the operators as 
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well. Related constructions have been considered by [SA] in the context of 
creation of spectral gaps. 
(2) For those who prefer tiling examples we now indicate how to view the con- 
struction above in this framework. Take a tiling dynamical system (see 
[LS2, Sol2| ) and replace one given tile T by a suitable homeomorphic image 
of T b indicated in Fi gure |[ We also indicated the next neighbor relations, 
showing that the resulting graph is just Gfi„ above. 



FIGURE 2. The tiling of Gfi 



Proof of Theorem uj. We first show that the condition is sufficient. Let u be 
an eigenfunction of A UQ associated to an eigenvalue E with suppw C B(x,r) and 
x G ujq. Then for any uj G fl every copy of P = B(x,r) A luq in Q A uj adds a 
dimension to the eigenspace of A U \Q belonging to the eigenvalue E. Let jpQ A w 
be the maximal number of disjoint copies of P in QAuj. Note that = : C 

is a upper bound for the number of points (and therefore the maximal number of 
copies of P) in B(0, 3r) n uj. 

This gives by a direct combinatorial argument that 

\ P Q Auj > -IpQ Alo. 

Thus, for arbitrary e > 

tr(X(-oo,g-e)(^|Q) ^(x^o^E+^jA^Q) _ 1 ftpUJ AQ 

\Q\ ~ \Q\ c \Q\ ■ 

Setting Q — with from a van Hove sequence and letting k tend to infinity, 
we get that p A (E - e) < p A {E + e) - 4p. Ase>0 

is arbitrary and v(P) > the 

desired implication follows. 
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Next we show the converse implication. Let E be a point of discontinuity of the 
function E i— > pa{{— oo, E}) and (Qk) an arbitrary van Hove sequence. We consider 



the distribution function j^n(A u , Q) of the measure Pq k - Proposition 3.4 shows, 
that jri-rn(A u , Qk) converges w.r.t. the supremum norm to the function E 
p A ((— oo, E]). Thus, for large k the jump at E of the function j^n(A UJ , Qk)(E) 
does not become small. More precisely we get 

dim (ker (A u \ Qk - £)) = lim(n(A,, Q k )(E + e) - n{A u , Q k )(E - e)) > c\Q k \ 

for a c > and all k £ N. Now let d 2 r A Qk = d 2rA Qk n Qfe denote the inner 
boundary of range 2r J 4 of Qfc. For a van Hove sequence (<5fe) we have 

dim £ 2 {d 2rA Q k ) = Hiel^ie d 2rA Q k n c^} 

< l^2r A +r(^)Qfc| 

- |B(0,r(w))| 

= £fc ' \B(0,r(u,))\ ' IQfc| 
with a suitable e k which tends to for k — > oo. For fc large enough we get that 
c > (k ' B(o,r(^))| • Thus, for large k the inequality 



dim (ker (A w | Qfc - £)J > dim e 2 (d 2rA Qk)- 
holds. Now let Wk the projection onto the inner boundary of range 2r^ of Qk- 



Then Propositon 3.5 shows that there exists an eigenfunction / of A^ such that 



Wkf = for k large enough. □ 

Remark 4.2. Let the conditions be as above. Then E is an infinitely degenerate 
eigenvalue of A^ for every The integrated density of states has a jump at 

E whose height is at least C~ 1 v{P). 

Acknowledgements: The authors would like to thank Uwe Grimm for helpful 
comments on the physics literature. 
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